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Secular instability 
in the spatial three-body problem 
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Abstract. Consider the spatial three-body problem, in the regime where one body 
revolves far away around the other two, in space, the masses of the bodies being arbi¬ 
trary but hxed; in this regime, there are no resonances in mean motions. The so-called 
secular dynamics governs the slow evolution of the Keplerian ellipses. We show that 
it contains a horseshoe and all the chaotic dynamics which goes along with it, corre¬ 
sponding to motions along which the eccentricity of the inner ellipse undergoes large, 
random excursions. The proof goes through the surprisingly explicit computation of the 
homoclinic solution of the first order secular system, its complex singularities and the 
Melnikov potential. 

1 Introduction 

The question of the stability of the Solar System is the oldest open problem in Dynamical 
Systems. A number of works have shown striking instability mechanisms in the three- 
body problems, e.g. oscillatory orbits and close to parabolic motion [Ale691 ISitbDl ILSSOal 
ILSSObl IGMS13i IDKdlRSli] , chaotic dynamics near double or triple collisions [BolObl 
IMoe89] (see also [M osni| l. But scarce mathematical mechanisms have been described 
regarding more astronomical regimes, which would be plausible for subsystems of solar 
or extra-solar systems. One of them [FGKRl^ shows the existence of global instabilities 
along mean motion (i.e. Keplerian) resonances. 

In this paper, we focus on secular resonances. Numerical evidence has long been 
suggesting that such resonances are a major source of chaos in the Solar system [LR931 
ILasOSl [FS89] . For example, astronomers have established that Mercury’s eccentricity is 
chaotic and can increase so much that collisions with Venus or the Sun become possible, 
as a result from an intricate network of secular resonances |BLF12] . On the other hand, 
that Uranus’s obliquity (97°) is essentially stable, is explained, to a large extent, by the 
absence of any low-order secular resonance [BLIP! ILR93] . It is the goal of this paper 
to provide a simple instability mechanism in the secular dynamics of the three-body 
problem. 
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Wniversitat Politecnica de Catalunya (Departament de Matematica Aplicada I), Spain, 
marcel.guardiaOupc.edu 


1 






























2 The secular Hamiltonian in the lunar regime 


Consider three point masses undergoing the Newtonian attraction in space. Call mo, mi 
and m 2 the masses and {qj,Pj)j=i,2,3 £ x the Jacobi coordinates. After having 
symplectically reduced the system by the symmetry of translations, the Hamiltonian 
depends only on {qj,Pj)j=i ,2 and equals 


J^Kep + Fi 


per, 


^Kep - 2^1=1,2 2p, W 


II 

2fj,j 

T? _ ( mam 2 

P"" “ Ik 2 || Ik 2 +mgi|| 


mi m 2 


where the reduced masses are defined by 

Mo = mo. Ml = mo + mi, + m]”^, ^2 ^ + m^ 


( 1 ) 


and coefficients cjj are the barycentric weights of the inner bodies: cjj = rrij/Mi (j = 
0,1). HKep is the integrable Hamiltonian of two uncoupled Kepler problems and induces 
a “Keplerian” action of the 2-torus, while Hper is the so-called perturbing function. 

In this paper, we consider the asymptotic regime (sometimes called lunar or hierar¬ 
chical or stellar) where the masses are fixed, while body 2 revolves far away around the 
other two. In particular, each of the two terms of KKep is negative, the outer eccentricity 
is bounded away from 1, and the semi major axes ai of qi and 02 of q 2 satisfy 


oi = 0 ( 1 ) <C 02 —00. 


The expansion of Fper in the powers of Ijf?! 11/11(7211 ^ 1 is 


F — — 

^ per — 


Pim-2 

Ik2|' 


(cos C) 


n>2 


|gl| 
1921 


(2) 


where the P^’s are the Legendre polynomials, and an = n > 2. Note 

that the perturbing function is of the order of l/o^. 

Standard normal form theor}0 shows that, for every integer k, there is a local change 
of coordinates o^^-close to the identity which transforms the Hamiltonian into 

P = PKep + Psec + 0(a2'^-'), (3) 

where Pgec is the secular Hamiltoniat^ of order k 

Psec = [ Pper Ml M 2 + O , (4) 

7t2 V«2/ 

obtained by averaging out the mean anomalies Ij at successive orders; is the orbit of 
the Keplerian action defined by PKep, parameterized by the mean anomalies ii and (.2 
of the two planets. The secular Hamiltonian descends to the quotient by the Keplerian 
action of , and induces a Hamiltonian system on the space of pairs of Keplerian ellipses 
with fixed semi major axes, describing the slow evolution of the Keplerian ellipses due 
to the perturbation. 

^See, in this context, [F^: mean motions do not interfere since they are not of the same order. 

^Secular, means century in Latin. This term governs the slow dynamics (as opposed to the fast, 
Keperian dynamics), on a long time scale, symbolically one century. 
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In this paper, we primarily consider the principal part of the Hamiltonian F of 
equation ([SD, i.e. F^ep + -^sec- 

We establish the phase portrait of the secular system, in particular with a (well 
known) hyperbolic periodic orbit located at inclined and nearly circular ellipses. The 
main results of the paper are the following. We prove that this periodic orbit posesses 
transversal homoclinic points (Corollary 14.2D and that the secular system posesses a 
horseshoe (Theorem 14.3p . and therefore chaotic motions. As a by-product, we prove 
that the secular system is non-integrable. The results we obtain are valid for any value 
of the masses of the three bodies. 

To prove these results, we take advantage of the fact that the first order in 1/02 of 
the secular Hamiltonian is integrable. We explicitely determine the homoclinic solution 
to the hyperbolic periodic orbit of this first order fLemma 13.ip . Then, we show that 
the Melnikov potential associated with this homoclinic orbit has non-degenerate critical 
points fProposition 14.1|) . from what it follows that the secular homoclinic solution splits 
for the full secular system, as well as for the full initial system. 

The result we present in this paper is a step towards proving the existence of unstable 
orbits in the three body problem in the lunar regime, for any value of the masses of the 
bodies, in the sense of Arnold diffusion. By “unstable orbits” here we mean orbits such 
that one of the three bodies undergoes a large change in the semi major axis of the 
associated osculating ellipse. This is explained more precisely in Section 14.11 


2.1 The quadrupolar and octupolar Hamiltonians 


The first two terms of the e^ansion ([2]) of Tper, after averaging, yield the quadrupolar 
and octupolar HamiltoniansO 


Tsec — /^1^2 (Tquad T (<^0 f^l)Toct) T O ^^5^ ’ 


(5) 


with 


and 


Tquad — 


ai 


8 02 ^ ( 1 - 622 ) 3/2 


((l5ei2 cos2g'i — 12ei2 — 3) sin2z -|- 3ei2 -|- 2) (6) 


^ 15 af 6162 

-^oct - - A ^ . r,. ^ /r. ^ 


64 a|(l-e2)5/2 


(7) 


cos gi cos g 2 


G! 


(5 sin2 i(—7 cos2 <71 -|- 6 ) — 3) 
— 35 sin2 gi sin2 i + 7 


sm gi sm g 2 cos ^ 


G' 


j (5 sin2 i{7 cos2 gi — A) + 3) 

-|- 35 sin2 gi sin2 i — 7 

this is a standard computation (see |Fej02| for the reduction to integrating trigonometric 
polynomials), and we have used the following notations: 


/ ... 

6j = eccentricities 

< gj = arguments of pericenters 

i = mutual inclination. 

\ 

®More generally, the n-th order term is called 2’*-polar, because, in electrostatic, this term is the 
dominating term of the potential of a system of 2" well chosen charged particles. 
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The Hamiltonians Tquad and Foct do not depend on the order k of averaging in the 
expression ([3]) (i.e., due to the special dependence of TKep and Fper in L 2 , these first two 
terms of the normal form are not modified by second and higher order averaging). 

Recall that we want to analyze the secular Hamiltonian F^ec for any value of the 
masses, so that parameters m*, /r* and Mj are just given constants satisfying ©• 

Let {£j,Lj,gj, Gj,6j, 2 denote the Delaunay variables. Jacobi’s classical elim¬ 

ination of the node consists in considering a codimension-3 submanifold of fixed, vertical 
angular momentum, and quotienting by horizontal rotations. The reduced manifold has 
dimension 8 , on which the Keplerian and eccentric variables {ij, Lj, gj,Gj)j=i ^2 induce 
symplectic coordinates. After averaging out the mean motions we are left with the 
symplectic coordinates {gj, Gj)j=i^ 2 , and the variables Lj may be treated as parameters. 

The variables Lj are given by Lj = Therefore, in the lunar regime we have 

Li ~ 1 and L 2 3> 1. Recall that the eccentricity of the ellipse is given in terms of the 
Delaunay coordinates by 



Since we want the outer body to describe non-degenerate ellipses, we even assume G 2 ~ 
L 2 . Since we are doing a perturbative analysis in L^^, we define the new variable 

r = C — G 2 , where G = 6 L 2 and J > 0 is a fixed constant. 


The coordinates 

(5I) <^1,7, T) = {gi,Gi, -g 2 ,G - G2), ( 9 ) 

are symplectic; we also call them Delaunay coordinates (as opposed to some radically 
different coordinates used later). Now these variables are bounded as 02 , L 2 —>■ 00 (recall 
that C ~ G 2 ^ 00 ). At the first order in L^^, the mutual inclination i satisfies 


2\2 


2. _{C-Gi- Gif _ (G^ -Gi-{G- T)") 

COS % — 




4Gf(C-r)2 


r2 


) • ( 10 ) 


Now we express the secular Hamiltonian in these new variables and expand it in inverse 
powers of L 2 . Recall that variables Lj are parameters, as well as the norm G of the 
angular momentum. 


Lemma 2.1. The secular Hamiltonian has the form 


with 


Fsec — oq + aiL2 ® (^Hq + L2 ^Hi + (do — 0 'i)L2 ^H2 + 0 {L 2 ^)^. 



2-5 



sin^ gi 


£! 

Ll 


( 11 ) 

( 12 ) 


In the coordinates (s'!, Gi, 7 , T), Lquad and hence Hq and Hi do not depend on 7 . 
We do not compute Hi explicitly, here. The reason is that Hi does not break the 
integrability of Hq and therefore does not play any role in the Melnikov analysis. In 
contrast, the Hamiltonian H 2 does break integrability (as it will follow from our study), 
and is computed in Section [5l 


Proof of Lemma \2.1[ Formula (l 8 |) implies 


J_ _ 
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Hence, using ([TO]) 


sin2 , = 1 - ^ + O (LJ‘) = 1 - - ^ + O (LJ‘) ^ (13) 

Besides, 

is a first integral of the first order approximation of -Fquad- Hence, expanding in powers 
of ^2 i^quad Can be written as 

i^quad(5i,Gi,r) = aiL2''(Fo(5i,G'i,r) + 2) + 0(L2-^) 

with 

_ LfM2lJ,^m2 

“ "8M2/i3(l_e2)3/2 

and 


Hq = (Se^ cos^ gi — Aef — l) sin^ i + 


= ei 


= e? 


2 (5cos2 5i- 4) ( +;^ + l 


GlJ ' Gj 


r 2 


2-5(l-^)sin^ffi 


£! 

Lr 


(using dH)) 

r 2 

— —2 (using (fT 0 |) and (fT3]) i 


Factorization m follows. 


□ 


3 The quadrupolar phase portrait 


According to ([ 6 ]), the quadrupolar Hamiltonian Fquad and thus Hq do not to depend 
on 7 . This happy coincidence (which does not repeat itself for the next order term 
Foct) makes Fquad integrable. This has been extensively used (see |Zhal4| and references 
therein). Here, we may thus think of T as a parameter. 

Complex singularities of solutions of Fquad are hard to determine in general. In 
our regime, the first order of the quadrupolar Hamiltonian, Hq, can be factorized as 
described in equation (fl^ (up to the additive constant —T^/L^), which dramatically 
simplifies our study. 

The Hamiltonian Hq is analytic on a neighborhood of the cylinder 


(gi,Gi) gT1x]0,Li[ 


in X M. Since it is vr-periodic with respect to gi, we may focus on the domain 
0 < <71 < TT, 0 < Gi < Ti, keeping in mind that the Delaunay coordinates blow up 
circular ellipses {Gi = Li). 

Hamilton’s vector field is 




Gi = 5 1 - 


5 1- 


Gi 


r 2 

G? 

1 - 


sin"^ <71—2 


r 2 

G? 


- 10 1 - 


Gl 


Gf 


sin^ gi 


sin 2 gi. 


The second component Gi vanishes if and only if (assuming Gi > 0) 


GiG{|r|,Li} or 51=0 (mod 7 r/ 2 ). 
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Gi separatrix 

Li 

r 

-^ 9i 

0 7r/2 tt 

Figure 1; Phase portrait of the reduced quadrupolar dynamics 
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We focus on the first two hyperbolic singularities, with Gi = Li, henceforth assuming 
that condition m is satisfied, with, say. 


r > 0 . 


(17) 


Lifted to the full secular phase space, these critical points become the periodic orbits 

^min,max(*>70) = > 7° + 7^ (*), , (18) 

where 7 ® is just the initial condition for the variable 7 and 

iHt) = -2^t (19) 

The whole circle Gi = Li corresponds to circular motion. Along this circle the Delaunay 
variables are singular. Thus, in a neighborhood of circular motion it is more reliable to 
use the Poincare variables 


i = y^ 2 (Li - Gi) cos gi 
rj = -y/2{Li - Gi)sing'i. 


( 20 ) 


Those secular coordinates are symplectic (a well known fact). The Poincare variables 
transform the Hamiltonian Hq into 


r2 1 

^o(C,7) = -^ + ^ 


2e^-(3-5^)r/" 


+ O 2 + 7^) 


( 21 ) 


and the line segment {Gi = Li} blows down to a single hyperbolic periodic orbit 

(C,7,7,r) = (0,0,7° +7^ (t),r) 


in the secular space. 

As it has occured above, in the (i^i, Gi)-plane there are two hyperbolic fixed points 
for gi E [ 0 , 7 r] and two more for gi E [vr, 27r]. On the line Gi = Ti, there are heteroclinic 
connections between them. All these objects blow down to the critical point (^, 7 ) = 
(0, 0) in Poincare variables. Moreover, there are two separatrix connections in the region 
Gi < Ti, one between the critical points with gi E [0, tt] and the other one between the 
critical points with gi E [tt, 27r]. Even if Hq is well defined for Gi > Ti, the corresponding 
solutions are spurious. In the secular space, we obtain a hyperbolic critical point at 
(^, 7 ) = ( 0 , 0 ) with two homoclinic orbits forming a hgure eight. 

The main technical goal of this paper is to show that these separatrices split when 
one considers the complete secular Hamiltonian of arbitrary order k. We focus on the 
homoclinic on which gi E [0, tt]. We use both Delaunay and Poincare variables, but need 
to keep in mind that the Delaunay variables are not defined on the secular space proper 
along circular inner ellipses (Gi = Li). Note the same orbit is homoclinic for the secular 
Hamiltonian m and heteroclinic after blow up (see (fT^ l. 

The expression of the energy (1121) yields the following characterization of this orbit, 
where the interval (s'™™,]0, 7 r[ is defined by the inequality (recall (|14l) l 

92 

sim^i > (22) 


and 



and 



(23) 
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Lemma 3.1. Assume conditions (USD -([HI). The Hamiltonian Hq given in (USD has a 
heteroclinic solution which tends to the periodic orbits and in (HHD in backward 
and forward times respectively. Its orbit is defined by the equation 


1-^J.in 9.=-. 


with range gi G ]0, 7 r[. Its time parameterization is given by 


with 


cos 51 (i) = - 


Giif) = |r| 


sinh A 2 t 


\Jx^ + (1 +X^)smh2 


3 ifi 

1 + sinh^ A 2 t (cosh ^20 


-1 


7 ^(i) = 7 ° + 7^(i) + 7 ^(i)) l^{t) = arctan (x ^tanhA 2 t) 


(24) 

(25) 

(26) 

(27) 

(28) 


where 7 *^ G T zs the arbitrary value of the 7 coordinate at initial time, 7 ^ has been 
introduced in (USD. 

In the 7 -component of the separatrix, the angle — 7 ^ is the (arbitrary) argument 
of the outer pericenter at the symmetry center of the separatrix 51 = 7r/2 (recall the 
change of coordinates ([SD). The term 7 ^(t) is the rotational part and the term 7 ^(t) is 
the transient part of 7 ^(t). 

Remark 3.2. The assumed condition (1161) implies that BtHq 7 ^ 0 in a neighborhood 
of the heteroclinic orbit. Therefore, one can rephrase Lemma l3.II as the existence of a 
heteroclinic orbit in each energy level HQ{gi,Gi,^,T) = h for any h < 0 (recall that by 
the expression of Hq in (I12p the heteroclinic orbits always have negative energy). 

Proof. On the separatrix, we have 

51 = lOT^-sm5f. 


Using 


and 


G\ 

, one can eliminate Gi by 

|r|\/5sin5i 


Gi = 


\/3-5 


cos^ 5i 


(29) 


and get a closed differential equation 

(1 “ 1(1 + X^)cos2 5i) ^1 - I cos2 51 


9 \ = 

where x has been defined in 


sm5i 


(30) 


and 


7li = 30V3r2(l-^^) >0. 










Note that, due to F 7 ^ 0 and ([ 22 ]), 51 / 0 along the separatrix. So, one can separate 
variables and, choosing gi = 7r/2att = 0 yields 


Alt 



_ singi dgi _ 

|(1 + x2)cos2 5(i) yr^^ycos^ 


The variable cos <71 G 
be dehned by 


Then 


Vs Vs 

V^V+X^) ’ V5(l+x2) 


is a coordinate on the separatrix. Let x 


cos gi = cos X G 


1 


Alt 



dx 

1 — (1 + x^) cos^ X 


1 



tan X — X 
tan X + X ’ 


with cos X 


y|cos5i. 


One can check that the constant A 2 defined in (|23p satisfies A 2 = y |AiX- The equation 
can be solved for tanx, which yields tanx = —xcoth( 742 t). This equality, jointly with 


cos gi = 



[3 ±1 

V 5 Vl + tan^ X ’ 


gives formula Using (12^1) . the Gi-component of the separatrix given in ([27)1 can be 
easily obtained. Finally, only the 7 -component remains to be computed. But, along the 
separatrix solution, 7 = drHo. Using ([T2I1 . (I2S)) and ([22)1 . 


7 =10 


4F 


(1 _ ^) si„7. - 2ij 

(1 - ^(1 + X^) cos^ 51 ^ - 2^ 


4F 1 F 

((1 + X^) cosh^ A2t - 1) Ll 


" lIx^aVE ^ (ESD follows. 


□ 


4 Splitting of separatrices 

Lemma l3.II shows that the Hamiltonian Hq in ([6|) has a heteroclinic connection (homo¬ 
clinic in the secular space for the Hamiltonian (I2ip i. 

The next term in the asymptotic expansion of the secular Hamiltonian given in 
Lemma O does not depend on 7 and therefore does not break integrability. For the 
amended Hamiltonian, according to classical perturbation theory, the hyperbolic periodic 
orbit and the homoclinic orbit persist. Moreover, the periodic orbit remains at (^, 7 ) = 
(0,0). Thus, in Delaunay coordinates, there are two hyperbolic periodic orbits, which 
are of the form 

f^min(^>7°) = {gmm + 0{e),Li,x° + x^{t) -FO(e),F) 

l^max(i>7°) = (5'max + 0(e),Fi,7° + 7ni) + G(e),F), 
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where £ = L 2 We can choose the time origin so that 


ye 

min,max 


(0,7°) = ( 5 . 


mm,max 


+ 0(e),Li,7°,r). 


The second order H 2 in Lemma [ 2.11 depends on 7 and therefore may (and will) break 
integrability. If we consider this Hamiltonian in Poincare coordinates (j20p . it posesses a 
hyperbolic critical orbit 0 (L 2 ^)-close to (^,??) = ( 0 , 0 ). 

To show that H 2 makes the stable and unstable invariant manifolds of this periodic 
orbit split, we use the classical Poincare-Melnikov Theory |Mel63] . More precisely, here 
we define the Poincare-Melnikov potential 




r+00 

(/)= / H 2 oZ^t,j^)- 

Jo 

+ [ H 2 oZ^t,j^) 


H 20 Z^^{t,-f^)dt 
-iL2oZTn(t,7°)df, 


similarly to [DG00] i. but taking into account the fact that the first order perturbation 
does not break integrability. 

In other words, we consider as an integrable system T-Lq = Hq + sHi and as a pertur¬ 
bation Til = H 2 (see Lemma l2.1l and recall that we have taken e = In Lemma l3.II 

we have obtained the parameterization Z^ of the separatrix (see (I25p i. For the Hamil¬ 
tonian "Ho, we have a parameterization of the separatrix Z^(t,7*^) = -|- 0(e). 

The time origin can be chosen so that 


Z^(0,7°) 



+ 0(e),7°,r 


(recall that by Lemma [5Tl Z°(0,7°) = (7r/2,Ty^sTs,7°,T)). 

In order to determine the splitting of separatrices, we consider the section gi = 1^12 
which is transversal to the flow locally in the neigborhood ot the unperturbed separatrix. 
It is still transversal to the perturbed stable and unstable invariant manifolds. We 
measure the splitting in this section. The distance between of the invariant manifolds is 
parameterized by 7 *^, the value of the 7 coordinate when t = 0 . 

Melnikov Theory ensures that the transversal homoclinic points in the section gi = 
TT!2 are = Lg^-close to the non-degenerate critical points of the Poincare-Melnikov 
potential. Expanding the potential in powers of e, we get £^( 7 °) = £( 7 *^) -|- 0{s) with 



H 2 oz^{t,y) 




r H 2 oZ^{t,j^) 

J —00 


H 2 o dt 

-/?2^°in(£7°)dL 


Moreover, since H 2 has a factor ei (see 0 ) it vanishes over the periodic orbits 
Thus, using (fT 8 ]l and ([ 2 ^ . the potential reads 



H2{g^,{t).G\{t),^\tlT)dt 
H2{g\{t),G\{t),^^ + ^\t)+^\t),T)dt. 


(31) 
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Melnikov theory then implies that, e-close to non-degenerate critical points of /I, there 
exist transversal homoclinic points in the section gi = 7r/2. The potential C naturally 
depends on parameters Li and T, which vary in the (non-empty) open set 

O = |(Li,r) : Li > 0,0 < r < LivW^} C ( 32 ) 

Let us write 

The next proposition shows that C necessarily has non-degenerate critical points, except 
maybe for exceptional values of the parameters. 

Proposition 4.1. There exists a non eonstant real analytic function : O —>• M such 
that the potential dal) is of the form 

+Z+i„re-*T'°. 

In particular, outside an analytic subset of O of empty interior, £+ does not vanish and 
thus C (as a function of'j^) has non-degenerate critical points. 

This proposition is proven in Section [5l where the function C'^ is computed explicitly 
(see formula ([38])). It is easier to describe the dynamics when the dimension is as low as 
possible, so let us carry out the symplectic reduction of the flow by time. Since 7 / 0, 
define the Poincare section = {7 = 7o} within some fixed energy level, and the 
corresponding return map, induced by the Hamiltonian (HID . 

V-fQ : y 

for which the Poincare coordinates (^,7?) may be used (see (l20p L This map has a 
hyperbolic fixed point Lf'^-close to the origin with one dimensional stable and unstable 
invariant manifolds. The classical Melnikov theory applied to the Melnikov potential 
obtained in Proposition 14.11 entails the following corollary and theorem. Note that the 
circle gi = 7 r /2 locally corresponds to the line ^ = 0 (see (l20l) f . 

Corollary 4.2. Fix (L5,r°) G O (defined in (1321) ) such that ^^0 po / 0 7 q 

be a non-degenerate critical point o/T^opp- Por L 2 large enough, there exists some 

7 q = 7 q -\-0{Lf^) such that the Poincare map V^* has a transversal homoclinic point in 
the line {.^ = 0} C . 

Recall that all Poincare maps are conjugate. Therefore, this corollary implies that 
there are transversal homoclinic points for V^g for all 70 G T. 

The transversality of invariant manifolds given by the corollary a priori refers to the 
secular Hamiltonian obtained after one step of averaging. Nevertheless, the conclusion 
holds for the Poincare maps of the secular Hamiltonian of any finite order. Indeed, higher 
order averaging only modifies higher orders of the asymptotic expansion of the splitting, 
while transversality at the first order was entailed by the first order of the expansion, as 
given by the Poincare-Melnikov potential. Hence, by considering the analytic set defined 
by the condition T^p ^p 7^ 0 at all orders of averaging, one gets the main result, where 
we restrict to some fixed compact set 

K cO = {(Li,r) : Li > 0,0 < P < Tiy^j C 

Theorem 4.3. Fix any 70 G T. There is an analytic set K° of full measure in K such 
that, for parameters in K°, for L 2 large enough, the Poincare maps : S.yp i->- 

S.p,p of the secular systems of all orders posess a horseshoe. 
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4.1 Splitting of separatrices for the non-averaged Hamiltonian 

As explained in Remark 13.21 for the secular Hamiltonian we have a hyperbolic periodic 
orbit at each energy level within a compact interval of such levels. These periodic orbits 
form a normally hyperbolic invariant cylinder. Corollary 14.21 implies that the invariant 
manifolds of this cylinder intersect transversally. 

In Section [5] we consider Li and L 2 as constants. One can lift the dynamics to the 
extended phase space, with the additional Keplerian coordinates (£ 1 , Li, ^ 2 ) A 2 ). Then, 
the cylinder gets enlarged by four extra dimensions: (Li,L 2 )) which are just constants 
of motion, and (^i,f' 2 ); which are performing a rigid rotation. 

This extended system, after rescaling, is just the first order of ([3]), that is 


Aq — Tkep T Fsec- 

Call A the cylinder of this Hamiltonian. Now one would like to analyze Hamiltonian 
([3]), that is, the full (non averaged) three-body problem. The first step is to prove the 
persistence of the invariant cylinder when L 2 is large enough, using that the reminder 
is for some A: E N. Since we are in a singular perturbation regime, the classi¬ 

cal theory of persistence of normally hyperbolic invariant manifolds [Fen72] cannot be 
applied directly. However, now there are results which deal with singularly perturbed 
problems and which can be applied to the present setting |Yan09[ fCdlLTlS] . Note that 
we can do global averaging up to high order. Therefore, one does not need to face the 
problems in [BKZllj and the obtained cylinder can be as smooth as needed. Call A the 
perturbed cylinder. 

In this setting, one can expect Arnold diffusion i.e., orbits whose action components 
drift by a large amount, uniformly with respect to large L 2 ’s. Since 

• Gi is prescribed by A and the homoclinic channel, 

• G 2 (resp. Li) is constrained by the conservation of the angular momentum (resp. 
the energy), 

one would expect to obtain orbits with drift in L 2 , that is in the semi major axis of the 
outer body. This is a remarkable feature, since semi major axes are known to be very 
stable when two of the three masses are very small [Nie96j . 

To obtain unstable orbits one usually combines two types of dynamics. The “inner 
dynamics” is the dynamics of Hamiltonian ([3l) in restriction to the cylinder A. The “outer 
dynamics” and is the so-called scattering map [DdlLSOS] . obtained as the following limit. 
Assume that the invariant manifolds of A split transversally along a homoclinic chanel. 
Consider a homoclinic orbit in the chanel, which is asymptotic to the trajectory of some 
point E A as t —>■ —00 and to the trajectory of some other point x+ E A as f ^ -|-oo. 
Then, we say that the scattering map S maps x_ to x+. hat 5 be a map (as opposed to 
a more general correspondance) is proved in [DdlLSOS] under general hypotheses which 
are satisfied here. Note that this map depends on the chosen homoclinic chanel and 
therefore it may not be defined globally -usually it is multivaluated. Provided that 
these two maps do not have common invariant circles, by iterating them in a random 
order, one gets pseudo-orbits which have the wanted unstable behavior. A shadowing 
argument then permits to approximately realize these pseudo-orbits as real orbits of the 
system. 

Understanding both the inner and outer dynamics is certainly not easy for Hamil¬ 
tonian Q. Concerning the inner dynamics, Jefferys and Moser |JM66j have used KAM 
theory to show that this cylinder contains quasiperiodic hyperbolic tori which form a 
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positive measure set. There should be a very rich dynamics in the complement of these 
tori in the cylinder. In particular, since the Hamiltonian restricted to the cylinder has 
three degrees freedom, there may exist Arnold diffusion in the cylinder itself. 

Concerning the outer dynamics, one needs first to prove that the invariant manifolds 
of the cylinder A split transversally and then derive some precise dynamical behavior 
for the scattering map, as in [DdlLSOS] . The analysis in the present paper allows us to 
perform only the first step. 

Theorem 4.4. Fix Lf > > 0. Consider the Hamiltonian dSI) with Li E [L-^ ,Lf] 

and L 2 > L 2 . If is large enough, 

1. there exists a normally hyperbolic invariant cylinder A, 

2. the invariant manifolds of the cylinder A intersect transversally along a homoclinic 
chanel, which is diffeomorphic to A, 

3. and there exists a scattering map S : A —>■ A associated to this homoclinic chanel. 


Note that, in the above statement, the cylinder is invariant in the sense that the vector 
field is tangent to the cylinder, but orbits may escape from its boundary (sometimes one 
rather refers to these manifolds as weakly invariant). 

Proof. As explained above, the persistence of the invariant cylinder can be obtained 
by the available results of persistence of normally hyperbolic invariant cylinders in the 
singular perturbation setting [Yann91 IGdlLTl,^ . 

For the other statements, we may use the results in Proposition 14.11 and [DdlLSO^ . 
Using the expression of Tgec in ©j one can split the Hamiltonian ([2D as T = Hq + Hi 
with Hq = F^ep — Fi'^ 2 FqvLa.d and Hi = F — Hq. Therefore, 

/a^ 

Hi = - /rim2(<To - fTi)Foct + O ( ^ 

\^2 

- //im2(<To - 0-i)Foct + o {Lf^) . 


and therefore satisfies Hi = 0 {Lf‘^). 

Let us abuse notation and assume that Hq and Hi are expressed in the variables given 
in d^D and that time has been suitably scaled. As we have explained the Hamiltonian Hq 
has as normally hyperbolic invariant cylinder. Recall that Li, L 2 and T are first integrals 
of Hq. Thus, the dynamics in this cylinder is foliated by three dimensional invariant tori. 
Let us define Yq = To(L ^ 2 ) -^ 2 ) 7*^) trajectory on the invariant tori L = L^, 

and r = in the cylinder with initial condition in the (^i,^ 2 ) 7 ) variables given by 


Since Hq is integrable, the stable and unstable invariant manifolds of the cylin¬ 
der agree, and form a homoclinic manifold. The latter is given (at the first order) by 
the homoclinic of the quadrupolar Hamiltonian obtained in Lemma 13.11 in the space 
( 5 i,Gi, 7 ,r). Recall that for Hq the variables Li and L 2 are first integrals. The 
dynamics of the variables £i and £2 is close to a rigid rotation and can be easily 
deduced. These homoclinic manifold can be parameterized by time t and the coor¬ 
dinates of the cylinder Y = Y{t, £ 1 , L^, £ 2 , L 2 ,'y^,T^). Note that Y can be asymp¬ 
totic to different points in the cylinder as t —?• ± 00 . There exist smooth functions 
{£f,£tn^) = (£f 7 ±)(^?,L?,£ 0 ,L 0 , 70 ,r 0 ) such that 


Y{t, £l Ll £l Ll 70, rO) - YQ{t, £f, L)’, £^, 7^, T^) 


zb rO /?di rO _,di ■p0\ 
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as t — > zboo. 

Then, to prove that the invariant manifolds split and that one can define the scat¬ 
tering map of the full Hamiltonian Hq + Hi, one may apply the results in [DdlLSOS] (see 
also [DdlLS06] i. Define the Poincare-Melnikov potential 

/ O r-\-oo 

HioY -HioYQ-dt+ / HioY-HioY+dt, 

-OO J 0 

where y stands for y(t, ^2)-^2’ 7 °; and yg^ for yo(t, L?, ^2) 7 ^; r°). 

Consider the function 


r I—)• C{f\ + WlT, T?, ^2 + ‘^2'T, ^2) 7^ + r'^) 


(33) 


where oj = (cji,a;2,wa) is the frequency vector associated with the torus L = L^, and 
r = Results in [DdlLSO^ imply that each non-degenerate critical point of this 
function gives rise to a transversal intersection of the invariant manifolds. The non¬ 
degeneracy of the critical point allows us to define a local scattering map. 

Using the formula of Hi, we see that 

£ = /: + c>(L2^) . 


where C is the Melnikov potential introduced in (1311) . Proposition 14.11 implies that, as 
long as £^o po 7 ^ 0 ; the function £( 7 ° — ojit) has non-degenerate critical points with 
respect to r. Thanks to the non-degeneracy, the function ()33l) has critical points, which 
are 0(£2 ^)"Close to those of £( 7 ° — ojit). Each critical point of (|33|) gives rise to a 
transversal intersections and to an associated scattering map. □ 

Theorem 14.41 implies that Hamiltonian ([HD fits in the classical framework of Arnold 
diffusion along single resonances: a normally hyperbolic invariant cylinder whose in¬ 
variant manifolds intersect transversally [DdlLSOOl IDdlLSOGl IBKZllj . Nevertheless, the 
results obtained in this paper are not enough to obtain unstable orbits drifting along 
the cylinder. Indeed, we cannot derive formulas for the scattering map and therefore we 
have no information about the outer dynamics beyond the fact that it is well defined. 

The reason is the following. To prove that the invariant manifolds of the cylinder 
split, it is enough to project them into a certain plane and see that they intersect 
transversally in this plane. But, in order to build unstable orbits, one needs more precise 
information on the scattering map. Namely, to derive a first order approximation of the 
scattering map, one needs to analyze how the invariant manifolds of the cylinder split 
in every direction (see |DdlLS08] i. in particular in the Li and £2 directions. This would 
give the size of the possible jumps that the scattering maps makes in these directions. 
That is, how far can be two points in the cylinder which are connected by a heteroclinic 
orbit. Since the mean anomalies ii and (.2 are rapidly oscillating in our regime, the 
transversality in the conjugate directions Li and £2 is exponentially small and therefore 
very difficult to detect. 

An intermediate step would be to tackle the 1-averaged Hamiltonian, where only 
the fastest mean anomaly ii has been averaged out. Then, one has a three degree of 
freedom Hamiltonian system with only one fast frequency. As a starter, we provide the 
expression of the analogue to the quadrupolar term: 


Ik2|| 



dil = 


14 
















/ ((4ei^ + l) sin^ 5 i + (l - ei^) cos^^i) cos^ i sin^ {v 2 + 92 )+ \ 

3//im2 10 e\^ cos gi sin gi cos i cos {v 2 + 92 ) sin (v 2 + 92 ) + 

Aa‘fr2^ ((1 - ei^) sin^^i + ( 4 ei^ + 1) COSTS'!) cos^ (^2 + 52) 

where V 2 = '^ 2 (^ 2 , L 2 ,-^ 2 ) is the eccentric anomaly. 

Analyzing how the invariant manifolds split either for the full three-body problem or 
for the 1-averaged Hamiltonian, and deriving formulas for the corresponding scattering 
maps, would be the major step towards proving Arnold diffusion in the three body 
problem in the lunar regime. 


5 Proof of Proposition 14.11 

In Lemma [ 3 .11 we have called {gi,Gi,'y^ = = T) the separatrix solution]^ 

We compute the potential defined in (I.Sip . To this end, we expand it in Fourier series in 
7 *^. Since H 2 is a trigonometric polynomial of degree 1 in 7 (or equivalently in 92 , see 
equation (UD), it can be written as 

H2 (51, Gi, 7, F) = Ht {gi , Gi, r)e'^ + (51, Gi, Fje"*^. 

Because 7 ^ = 7 *^+ 7 ^(t)+ 7 ^(t) depends linearly on 7 *^, too has only two harmonics: 

/:(7o) = /:+e*T'°+/:-e-'^°, (34) 


where 


= 


/ OO 

Hf {g\ (t ), G’l (t ), 

-00 


Since £ is a real function, it is determined by, say, the positive harmonic C~^ = C~. 

As a first step, we parameterize the separatrix by gi instead of t, using the following 
lemma (where we omit the upper index h). 

Lemma 5.1. The following identites are satisfied on the separatrix given by Lemma \S.l\ : 

( 



The proof of this lemma is a direct consequence of Lemma (3.II and formulas ([8]), (HOI) . 
m- We can use this lemma to express the function = H 2 on the separatrix as 
a function of gi. 

^The variable t used in section|3]is a rescaled time, since in Section|3]we have simplified the quadrupo- 
lar Hamiltonian by some multiplicative constant, when replacing Fquad by Hq. 
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Lemma 5.2. The function J-^ can be written, on the separatrix, as a funetion of gi as 
= \{Fi + iF 2 ) with 


1 -^ 


TT7 


\/l—i (l+X^) cosgi^ / _ 

-^1 COSgl 1-I^COSgf 


^ „ Ji-^{i+x^)cos^gi 

-^2 = ^2 - / =- 

^ 1-3 cos^ 31 


5 cos2 31(9-11 cos^ 31) I 

-^—5- 2 -r 

1—^ cos^ gi 


3lf 


and 


Cl = 


105 af 62 


32 a|(l-ei)5/2 

Proof. We have 


1 - 


5r 


2\ 3/2 


3LlJ 


^-5^- (14-11 


and C 2 = 


cos^ 51 


15 /^of 62 


64 V 3 a|(l - 62)5/2 


1 - 


5r2\ 1/2 

sZfJ 


H 2 = Aoct^i (cos 5 i cos 7 A + sin 51 sin 7 cos z B) 


where 


and 


Thus, 


A= jj {p sin2 i(—7 cos2 51 + 6 ) — 3) — 35 sin2 gi sin2 i + 7 
B = ^{b sin2 i{7 cos2 51 — 4) + 3) + 35 sin2 gi sin2 i — 7 


-Tort, — 


15 a? 


62 


64 a|(l-62)5/2- 


-(- _ ToctSl 


= 


{A cos 5 i cos 72 + -B sin 51 sin 72 cos i) 

+ i{A cos 5 i sin 72 — -B sin 51 cos 72 cos i) 


which we want to express in terms of 51 . Using Lemma l5.11 the functions A and B can 
be written as 


A = 


5 


31? TT; 


B = - 


I cos^ gi 

5 r2 


(9 — 11 cos2 gi) — 7 


I ^ 1—5"^“^— (5—11 cos2 5i) + 7. 

3 Lf 1-| cos2 31 ' ' 


Let 


Then, 


Bi = 6 i(cos 5 i COS 72 T + sin^i sin 72 coszB) 
S 2 = 61 (cos 5 i sin 72 T — sin 51 cos 72 cos iB). 


Si = - 


i 4 r 

LiX 


1 - 


5r2 


- |(1 +X^)cOs2 5ri 
1 - I cos 2 gi 


and 


B 2 =- 


VW T ^1- i(l + x2)cos2 g^ 


cos 5 i 1 - 


-1 11 r 2 

5 ^ “ TI? 

3 ~i 5 r 2 

'^1-31? 


cos 2 5i 


3 Lix 


1 - I cos2 gi 


5 r2 C0s2 5r]^(9 — 11 COs2 5l) 21 r2 


3L2 


1 - ^ cos 2 gi 


Ll 


r 2 


+ — - 5^ - 14 - 11^ cos 51 


Ll 


□ 
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Now we compute J-^ as a function of t. Recall that the constant A 2 was defined in 


Lemma 5.3. The functions J-i and J -2 on the heteroclinic (cf. Lemma 
written as functions of t = A 2 t as 


sinh? 


7+6sinh^ r 


1 +sinh^ T x^ + (l+X^) sinh^ r 

-^2 = 


5 + 


sinh^ ' 


X2+(1+X^)sinh"= T \L(x 


p 2 9x^ + (f+9x^)sinh^T 3 11 

cosh^ T 5 Zfy 


with 


Ci = - 


15 af 62 


16^10 a| (1 - e|)5/2 Li [ 3 Lf 


1 - 


5r2 


Proof. We compute each term as a function of r. By (|26D . we have that 

X 


- ^(1 +X^)cOS5r5’(t) = — 


and 


So, 


+ (1 + x^) sinh^ r 
5 0 2 1 + sinh^ r 


\/l - 3(1 +x“)cos 9 f(t)cos 95 '(() 


1 / 3 sinh T 


1 - I cos 5 ? (t) 


For the last term, we have 




X V 5 1 + sinh^ r 
1- |^)x^+Msinh2r 


1 - 


5 r2 
3l? 


(x^ + (1 + X^)sinh2 


where 




7L? 


Using the dehnition of x in 


we have that 


1 - 


5r2 

3Lf 


X = 


Therefore, 


1 - 


1 - - 


11 

“I? 


|!-cos2fi(?(f) I = 


3 1 5 r2 

3Zf 


2r2 

sZf 

4r2 

" = 7lf 


2r2 


7 + 6 sinh^ r 


21 L 2 1^1 _ I^ j X^ + (1 + X^) sinh r 


Putting all these formulas together and defining 


Cl = 


21x 


Cl 


3r2 




con be 


□ 


we complete the proof. 
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The formulas of Lemma 15.31 allow us to compute the Poincare - Melnikov potential 
We split it as as C~^ = Ti + i £2 with 


~ 2 




(35) 


We first compute 


II 

sinh(T2f) 

7 + 6sinh^(T2f) 

_oo 1 + sinh^(T2f) 

X^ + (1 + X^) sinh^(T2f) 

Cl 

sinh r 

7 + 6sinh^T 

A 2 . 

/_oo 1 + sinh^ r x^ 

+ (1 + x^) sinh^ r 


I -?T L 

^ dt 


using the residue theorem. First note that if we look at this integral along the path 
T = —in + u, (T G M, instead of the real line, we get 


• —27r+R 


sinh T 

1 + sinh^ r + (1 + X^) sinh^ r 


7 + 6sinh^T 

e ^ 2^1 dr 


= -e 

So, if we define the function 

fii'r) = 


/ 

Jr 


sinhfj 


7 + 6 sinm a 


1 + sinm a + {1 + x^) sinm a 


e da. 


sinhr 


7 + 6 sinh^ r 


1 + sinh^ r + (1 + X^) sinh^ r 
Cauchy’s integral formula shows 

27rz 




Ci = - 


1 + e ^ 2 + 


2r-7r 

AoLi 


■ Residues 


(36) 


(37) 


where the summation stands over the residues of the function fi in —n < Imr < 0; the 
negative sign comes from the index of the curve. 

Lemma 5.4. The function fi defined in (IMl) has three singularities in the region —n < 
Imr < 0 given by 

/ — • TT 

«0 = -*2 

or = —i arcsin 


af = —i { n + arcsin 




and the associated residues are 


R-es (/,ao ) =--^e 


Res (f,a7) = 

2 y/Tfix 


2r • V 

2 — - —7Y Rrcsin —^ 


g ^2-^1 x/l + X^ 




2r 

2 - T-PI 


A 2 Lf 


TT— arcsin 


V'l+x2 


Before proving the lemma, we proceed to deduce the value of the real part of the 
Melnikov integral. The lemma shows that the sum of formula (1371) is 


Residues 


2r - 


A 2 LI 


e ^ + 


7 + X^ 


2 ^1+ X2 


or • V 
—:—T arcsin , ^ ^ 
g ^2^-2 Vl+x2 


1-e 212^1 
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hence the function £i introduced in (f35]l satisfies 


A = - 


27ri 


Cl 


2 r^e ^ + 


1 + e 




Ao 


+ 


A 2 LI 
Lhx!p" 


OF 

0 arcsin ■ 


1 — e 


2r . 


Proof of Lemma \5.4\ The singularities are the solutions of 

1 + sinh^ r = 0 or either + (1 + X^) sinh^ r = 0. 


\ 


The first equation is just 1 + sinh^ r = cosh^ r = 0 and the only possible solution in 
—IT < Imr < 0 is Oq = —mj^. For the second equation, one can take r = ia. Then, it 
is equivalent to 


sin (T = ± 


X 




which gives the other two singularities. 

Now we compute the residues. We start by Oq . We compute the Laurent series of 
each term in /. From the fact that 


sinhr = —i (1 + 02 {j + ittIT)) 


one can easily deduce the following, 


1 

cosh^ T 
7 + 6 sinh^ r 
X^ + (1 + X^)sinh^T 


“7 , • I^\2 ^ 2 ( 1 " + iT^IT)) 

(T + ITT /1)‘^ 

1 + 02{t + iTT/2) 

-1 + 02{t + iir/T). 


Therefore 

sinh r 7 + 6 sinh^ r 

1 + sinh^ r + (1 + X^) sinh^ r 

For the exponential, we know have that 


i 

(t + iiT jT)'^ 


(1 + 02 {t + f7r/2)) . 




2]^ _ ttF 

i . 2 ^ + i7r/2) + 02 ( 1 " + fvr/2). 

A2L^ 


From these two last expansions, we obtain the residue for Og = 

Now we compute the other two residues. We compute them at the same time. Note 
that for z = 1, 2, 

sinh a” = —i — , ^ 


and 


Therefore 


cosh 


1 

7T+7’ 


cosh 02 = 


1 

Vl + x^ 


7+^2 

7 + 6 sinh^ r =-^ + OAt — iir/T) 

1+X2 

X^ + (1 + X^) sinh^ T = -2xz(r - oj") + C> 2 (r - oj") 
X^ + (1 + X^) sinh^ r = +2xz(r - o^) + 02(r - o^) 
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and we have also 


■ r.^ 


A 2 -L 


1 =e ^2^1 Vi+x2 + ) 


■ rA 


( TT—arcsin-^ 


+ C)i(t - ). 


g A 2 -L 1 _ g ^2^1 

With all these expansions, one can easily deduce the last two residues. 
Analogously, the function £2 introduced in (1351) . satisfies 

C 2 27ri 


C.2 = - 


Ao 


Residues, 


1 + e ^ 2^1 

where the sum stands over all residues in —tt < r < 0 of 
2r 


/2('r) =- 


— % - TT' 

O 242^1 


cosh T 


21 r2 


sinh^ r 


+ {l + X^) sinh^ r V L\x 


r2 9x2 + (f + 9x2) Sinh2r _ 3 _ 


cosh^ r 


l! 


□ 


Lemma 5.5. In the region —tt < Imr < 0, the function /2 has the same singularities 
as fi, given in Lemma\5.4\ Moreover, the associated residues are given by 


Res (/2,ao ) = 

-i ( 21 - 81 ! 

+ 24- 

Res {f 2 ,af) = 

,•3 X / 


'^5(1+x2)3/2 ( 

Res {f 2 ,af) = 

,•3 X / 


'5(l+x2)3/2 ( 


rf_ 

IfIfP 


e ^ 


.. , TT— arcsin , ^ ^ 

V vWW 


Hence, the function £2 introduced in ([35]l satisfies 

/ 

C 2 27ri 


C.2 = - 


A.2 


1 + e ^ 2 rf 


-^21-8^ + 24 




e ^^2 + 


p 2 \ “ -77 arcsin 

(l+x"2)3/2 - 7 I e 2 


,•3_X 


Vl+x2 _|_ 


V 


—I 


,•3_X 


llir-T e ^2 


My I TT—arcsin , 


/ 


^5(1+;^2)3/2 ^ 

Proof of Lemma 15.51 We use the expansions obtained in the proof of Lemma 15.41 For 
the last two residues, it is enough to use also that 

r2 9x2 + (f + 9x2)sinh2r 33 F^ ^ 

If?--= Ti2+oh-“.) 

for i = 1,2. Then, using the expansions of Lemma 15.41 we can deduce the two last 
residues. 

Now we compute the residue at r = oq. We split /2 into three parts f 2 = hi + h 2 + h^ 
with 


hi{T)= (¥-5^^ 


5 ^ I cosh 5 


I - TT-) 

A2Li 


<h 2 (r)= -1(14-11-^ 


sinh^ T 


— 't -?r' 

AnL? 

e 


L? y x^+(l+X^) sinh^ T coshr 
2 r 


h.(T)= 9xW(f+9x") sinh" ^ 2 e ^ 

3V ) L'fx^ X^+(l+X^) sinh2 r cosh^ r 
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Functions hi and /i 2 have a pole of order one at r = Oq so one can proceed as for the 
others singularies. However, hj, has a pole of order 3 and therefore, we need to compute 
the expansion up to order 3 of each term in h 2 - 
We use the expansions 


jsinhr =-i - + Oi{y) 

I cosh r =-iy - ly"^ + 05 {y) 

where y = r + We use this notation throughout the proof. 
We start with hi and h^. It can be easily seen that 


irr 


Res(/ii,ao) = ^ e 


Res(/i2,ao) = ^ 


' ML'i 


For /13 we use the following expasions 


9x^ + (t + 

+ (1 + x^) sinh^ r 


, 2 12 /12 33 ., , ., , 

smh^ r = —- - — + —X V + Oi{y) 


1 


cosh^ r 
»■ 2r _ 




V 5 ' 5 

y^ 


— 3 ' 1 + ^>4(2/) 


= e [l-i 


2T 


2T^ 


^A2LY aIl\ 

Putting together all these expansions one can deduce the residue 


y^ + 0^{y)]. 


T> -\ ■ 1^33 2 24 r2 \ - -r 


e ^ 2 "^! 


Now it only remains to add the three residues to obtain the residue of /2 at r = Oq . □ 

Gathering what precedes, and letting 


ttF 


a = 


F = 


X = 


A 2 LI^ * Li 

we obtain the following analytic expression; 


— , and p = — arcsmx, 


/ 


£+ = - 


27ri 


A 2 (1 + e-2“) 


Cl 


-^e-“+ 
7+X^ 


^p-2/3 fl _ p-2q 


+ 2 ^sin^e-^P ( 1 -e 

/ i 


) 


+C 2 


V 


~ c: 


21 _ SF^ + 24^ j e-“+ 

'llf 2 - 7 ) e-2/3+ 

\lf^ - 7 I e-2(“-/3) 


5 (l+x2)3/2 

1 -U —_^_ 

V ^5(1+x2)3/2 


(38) 




In order to check that £'*' is not consant, notice that, asymptotically when F tends to 0 
(and Li is kept constant), a = 0 (F“^) while /3 = 0 (F“^), so that the term in Ci 
dominates the others: for some C 3 (Li) and ( 74 (Li) independant of F, 


£+ ~ exp (-2C'4(Li)r + o(r- 2 )). 
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Since the dominant term of the right hand side is a non constant function of F, Propo¬ 
sition O is proved. 
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